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ABSTRACT 

We analyse the global structure of the old stellar discs in 34 edge-on spiral galaxies. The 
radial and vertical exponential scale parameters of the discs are obtained by applying 
an improved two-dimensional decomposition technique to our /-band photometry. We 
find a clear increase in the disc scaleheight with maximum rotational velocity, in 
accordance with observations of the stellar velocity dispersions in galaxy discs. The 
range and maximum of the intrinsic flattening of the disc light seem to increase with 
both maximum rotational velocity and total Hi mass. We use the disc flattening to 
estimate the disc contribution to the maximum rotational velocity, resulting in an 
average of 57±22 percent. The disc light distributions arc further investigated for the 
presence of radial truncations. We find that the radial light distributions of at least 20 
spirals are truncated, corresponding to 60 percent of the sample. For small scalelength 
spirals, which are the most numerous in the local Universe, the results suggest that 
the average ratio of disc truncation radius to disc scalelength is at least four. 

Key words: galaxies: spiral - galaxies: structure - galaxies: fundamental parameters 
- galaxies: photometry - galaxies: stellar content 



1 INTRODUCTION 

The study of the light and mass distribution of galactic discs 
is of fundamental importance for our understanding of the 
formation and evolution of spiral galaxies. It is well estab- 
lished that globally, the radial light distribution of galactic 
discs is exponential in nature (Freeman 1970), over a wide 
range of central surface brightness (de Jong 1996b; de Blok, 
van der Hulst, & Bothun 1995). This radial exponential dis- 
tribution can be understood in the framework of the collapse 
theory of disc galaxy formation if the stellar mass-to- light ra- 
tio is approximately constant with radius (Fall & Efstathiou 
1980; van der Kruit 1987; Dalcanton, Spergel, & Summers 
1997). Perpendicular to the galactic plane, the light distri- 
bution of the disc is also approximately exponential (van der 
Kruit & Searle 1981a,b, 1982, hereafter KS1-3; Wainscoat, 
Freeman, & Hyland 1989; Shaw & Gilmore 1990; de Grijs, 
Peletier, & van der Kruit 1997). While the radial structure 
of the stellar disc is probably determined during disc galaxy 
formation, the vertical structure largely reflects the secular 
disc evolution (Wielen 1977; Bahcall 1984; Lacey 1984). To- 
gether, the radial and vertical disc structure determine the 
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intrinsic disc flattening. Accordingly, disc flattening is gov- 
erned by both disc galaxy formation and evolution and may 
offer additional insight into the physics of these. 

The disc flattening can only be studied accurately in 
edge-on galaxies, where the radial and vertical distributions 
are essentially independent (KS1). Previous studies of large 
samples of edge-on disc galaxies pointed out that, on aver- 
age, the disc light distribution slowly flattens from earlier to 
later morphological type (Guthrie 1992; Kudryaet al. 1994). 
These studies adopted a statistical approach, using isopho- 
tal axial ratios. A similar trend is observed in a sample of 47 
edge-ons (de Grijs 1998), using the ratio of the exponential 
scalelength to the exponential scaleheight of the disc light 
(see also Schwarzkopf & Dettmar 2000). In this sample, a 
further comparison of the disc flattening with other global 
parameters, such as rotational velocity or H I content, may 
reveal additional trends concerning disc galaxy formation 
and evolution. Besides this, the flattening of the disc light 
can be used to estimate the disc contribution to the rota- 
tion curve (Bottema 1993). This allows one to study the 
dark matter content of spiral galaxies and provides an inde- 
pendent test of the so-called maximum disc hypothesis (van 
Albada et al. 1985; Freeman 1992). 

A related study of the stellar kinematics of galactic discs 
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by Bottema (1993) shows that more massive spirals tend to 
have larger stellar velocity dispersions. Together with the 
increase of the stellar velocity dispersion with age (Wielen 
1977), this suggests that more massive discs are more dy- 
namically evolved. If true, then more massive discs are also 
expected to have larger scaleheights. This can easily be ver- 
ified observationally in edge-on spiral galaxies. 

Another global parameter describing the radial light 
distribution of discs is its outer edge. This truncation of 
the disc light occurs at a radius fl max , which varies from 3 
to 5 scalelengths among disc galaxies (van der Kruit 2001). 
It is most easily found in edge-on spirals because of their 
higher surface brightness (KS1-3; Barteldrees & Dettmar 
1994; Pohlen et al. 2000b; Schwarzkopf & Dettmar 2000; 
de Grijs, Kregel, & Wesson 2001). In less inclined spirals, 
for which azimuthally averaged radial light profiles are often 
studied, the non-axisymmetric component (e.g. due to spiral 
structure) tends to smooth out a truncation present in the 
old disc. This effect was first noted by van der Kruit (1988); 
for 16 face-on spirals, of which 15 did not show any sign of a 
truncation in the radial light profile, he found that the three 
outermost isophotes were much more closely spaced than the 
inner ones, providing clear evidence for truncations. The ori- 
gin of the truncation is still unclear. If the truncation is also 
present in the disc mass distribution, it may be caused by the 
angular momentum distribution of the proto galaxy (van der 
Kruit 1987). This would imply that Hi observed beyond the 
stellar disc (e.g. Broeils & Rhee 1997) has been accreted (cf. 
Bottema 1996). Another possibility is that the truncation 
corresponds to the radius at which the gas density drops 
below a threshold density necessary for star formation (Fall 
& Efstathiou 1980; KS3; Seiden, Schulman, & Elmegreen 
1984; Kennicutt 1989). 

In a previous study, de Grijs (1998) analysed BVIJK' 
photometry of a statistically complete sample of nearby 
edge-on disc galaxies. He used profiles taken at suitably 
chosen positions to determine both the radial and verti- 
cal scale parameters (hereafter ID method). The determi- 
nation of the radial scalelengths was based on radial pro- 
files, taken at some distance above the plane to minimize 
extinction. Still, the uncertainty in this analysis is consid- 
erable. First, spiral galaxies may have a surface brightness 
distribution which shows significant local deviations from as- 
sumed smooth model distributions (Shaw & Gilmore 1990; 
de Jong 1996a). Since the ID method uses a small portion 
of the galaxy, it is susceptible to these local deviations. Sec- 
ondly, the ID method does not consider the presence of a 
radial truncation. Finally, an exponential was used as a fit- 
ting function, thereby neglecting the effect of line of sight 
projection. These omissions have probably caused a serious 
error in the disc scalelengths. 

Here, we re-analyse the /-band photometry of a sub- 
sample of de Grijs's sample, consisting of 34 spiral galaxies. 
In the first part of this paper, we use a two-dimensional 
(2D) least-squares algorithm to study the exponential scale 
parameters, and thereby the disc flattening. This 2D method 
minimizes the effects of a radial truncation and takes the line 
of sight projection into account. In addition, it allows for a 
bulge-disc decomposition and adopts a mask to minimize the 
effect of dust extinction. To assess the effect of dust extinc- 
tion outside this mask, the B-band data are subjected to 
the same analysis. The 2D method was briefly introduced 



by de Grijs et al. (2001), where the results of a detailed 
study of the radial truncations of four spiral galaxies were 
presented. In the second part of this paper we report on a 
similar study of the radial truncations for our sample in the 
/ band. Finally, the observed truncation radii are compared 
to those predicted by the scenarios proposed for the ori- 
gin of the truncation. This paper is organized as follows. In 
Sect. 2 we present the sample. The 2D least-squares method 
is discussed and tested in Sect. 3. In Sect. 4, the results of 
its application to the /-band photometry are presented, to- 
gether with the results of our study of the disc truncations. 
We further analyse and discuss the results in Sect. 5. Our 
main findings are summarized in Sect. 6. 



2 THE SAMPLE 

The parent sample (de Grijs 1998) contains 47 spirals and 
lenticulars and was selected from the ESO-LV catalogue 
(Lauberts & Valentijn 1989) according to the following cri- 
teria: 

(i) inclination, i > 87°; 

(ii) blue diameter, Df 5 > 2'. 2; 

(iii) Hubble type ranging from SO-Sd; 

(iv) non-interacting 

From this sample we selected all regular spiral galaxies and 
excluded those that are clearly lopsided or warped. Impor- 
tant properties of the 34 selected spirals are summarized 
in Table 1, columns (l)-(6) ' . We note that although the 
sub-sample covers a large range in maximum rotational ve- 
locities («max = 50-400 km s _1 ), it is dominated by spirals 
of intermediate- to late- type. 

We checked the completeness of the sub-sample with 
the V/Vmax test (Davies 1990). Based on the blue angu- 
lar diameters given by the ESO-LV catalogue, the average 
V I Fmax is 0.58 ± 0.06. This is slightly larger than the av- 
erage of 0.5 expected for a randomly distributed sample. In 
the remainder of this paper, one should keep in mind that 
our sample is of modest size and may not be statistically 
complete. 

While our J and if'-band data more closely resemble 
the underlying disc mass distribution than the optical BVI 
data (Rix & Rieke 1993), the quality of the former is poorer 
(de Grijs 1998) and will not be used. Of the optical pass- 
bands, the light in the / band is the least affected by dust 
extinction and the most dominated by the old stellar pop- 
ulation (e.g. Chariot & Bruzual 1991). We will therefore 
primarily use those data to study the global structure of 
stellar discs. The data reduction and sky subtraction have 
been discussed by de Grijs (1997, 1998). 



3 THE TWO DIMENSIONAL METHOD 
3.1 The model 

The spiral galaxy luminosity density is assumed to be com- 
posed of a disc and a bulge component. The disc is assumed 

t All distance dependent parameters in this paper are calculated 
using H =75 km s _1 Mpc -1 
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to be axisymmetric, transparent and, in first instance, with- 
out a radial truncation: 

L(R,z) = L e - R/hR e - z/h % (1) 

where (R, z) are the cylindrical coordinates, Lq is the central 
luminosity density and /ir and h z are the radial scalelength 
and the vertical scaleheight (in terms of the vertical scale 
parameter zo; h z — 0.5 zo). While other functionalities have 
been proposed for the vertical distribution, most notably the 
locally isothermal (KS1) and sech(z) distributions (van der 
Kruit 1988), the exact choice is unimportant for our pur- 
pose (Sect. 3.2). The vertical exponential behaviour is taken 
to be independent of radius following de Grijs & Peletier 
(1997). They showed that, in the discs of most of the spirals 
in our sample (at least all late-type galaxies, which is about 
two thirds), the scaleheight does not increase with radius. 
It is further assumed that the galaxies are perfectly edge- 
on. As KS1 demonstrate, the error of the scale parameters 
due to this assumption is negligible as long as i > 86° (see 
also de Grijs et al. 1997), which is the case for all the galax- 
ies in our sample. With these assumptions the disc surface 
brightness distribution is then given by (KS1): 

T, diBC (R',z) = S (R'/h R ) K^R'/hn) e - z/h \ (2) 

where R' is the projected radius along the major axis, Eo = 
2/irLo is the projected edge-on central surface brightness 
and K\ is the modified Bessel function of the first order. 

The surface brightness distribution of the bulge is mod- 
elled either by an exponential: 

T, hulge (R',z) = 5.360 S c e- 1 ' 679 r/r % (3) 
or by an r 1//4 law: 

S bulgc ( J R',z) = S ce - 7 - 669 ^ 1/4 - 1 ), (4) 

here r is the projected radius; r — \J R' 2 + (z/q) 2 , r c is the 
effective radius, E c the effective surface brightness, and q the 
bulge axial ratio (the ratio of the bulge minor to major axis 
diameter). The best-fitting bulge, in the reduced \ 2 sense, 
is adopted. To account for the effect of seeing, the model 
bulges are convolved with a two-dimensional Gaussian PSF. 
For the disc component the seeing effect can be neglected 
(Sect. 3.3). The central positions and position angles of the 
galaxies, determined by fitting ellipses to the /-band images 
(Sect. 4.1), are fixed parameters in the fit. Thus, the adopted 
model has at most six free parameters; So, hn, h z , E e , r e , 
and q. 



3.2 The approach 

Since the studies of van der Kruit & Searle (KS1-3), several 
techniques have been developed to model the global struc- 
ture of stellar discs in edge-on spiral galaxies (Kylafis & 
Bahcall 1987; Shaw & Gilmore 1989; Xilouris et al. 1997). 
The application of these techniques to the region near the 
galactic planes was found to be problematic, mainly because 
of three reasons; (1) the presence of large amounts of dust, 
(2) the presence of additional components, e.g. in the form 
of spiral structure or a young thin disc, and (3) the limited 
accuracy, of about 1-2 degrees, with which the inclination 
can be determined. 



Photometric studies of large nearby spirals which in- 
clude a radiative transfer description, have pointed out that 
dust is more confined to the plane than the disc light; 
/i z ,dust ~ 0.5 /i z , stars (Kylafis & Bahcall 1987; Xilouris et al. 
1997, 1999). Considering the relatively small velocity dis- 
persions of young stars (e.g. Wielen 1977), the same is true 
for young stellar populations. Hence, if we assume that the 
range in z/h z affected by dust extinction and young popula- 
tions is similar among galactic discs, it becomes possible to 
mask out that region in a systematic way. Indeed, the verti- 
cal / — K colour profiles of the galaxies in our sample, which 
trace the amount of dust extinction, show that in all cases 
the region \z\ > 1.5 h z is only marginally affected by dust 
extinction (de Grijs et al. 1997). Therefore, to solely study 
the old stellar population, we choose to mask out the region 
corresponding to \z\ < 1.5 h z . Since the differences among 
the various functionalities proposed for the vertical light dis- 
tribution are small outside this mask (van der Kruit 1988, 
his fig. 1), the choice for the vertical distribution becomes 
unimportant. Still, we will further assess the presence of dust 
extinction outside the mask by (1) testing the 2D method 
on 'dusty' artificial images (Sect. 3.3) and (2) repeating the 
2D fits in the B band to make a comparison with the /-band 
results (Sect. 4.1). 

The 2D decomposition uses the non-linear least-squares 
technique developed by Marquardt (1963). The fit of 
Eqs. (2-4) is done in the linear regime, to the masked version 
of the sky-subtracted image (the data are not binned as in 
many previous methods, e.g. KS1-3, Shaw & Gilmore 1989). 
To ensure that the fitting algorithm minimizes the relative 
difference between our model and the data, each pixel is 
weighted inversely proportional to the value of the model 
surface brightness distribution of the disc at that pixel (cf. 
de Jong 1996a). As a consequence of this weighting, regions 
of the lowest signal-to-noise carry a large weight in the fit. To 
prevent this side effect, the data outside the region defined 
by the la contour of the model image are excluded from 
the fit (the la level is the standard deviation of the back- 
ground) . By using this criterion on the model instead of the 
data, there is no discrimination in favour of positive noise 
peaks. Since the adopted weighting scheme and the fitted re- 
gion depend on the initial estimate of the model, an iterative 
application of the algorithm is required (see Sect. 3.3). 

The light distribution of galactic discs may have a grad- 
ual or sharp radial steepening in the outer parts, indicating 
a possible truncation (e.g. KS1-3). If present, such a region 
will be excluded from the 2D fit by using an outer fitting 
boundary. This boundary will be assigned by eye, at a pro- 
jected radius within the radius at which the radial light pro- 
file starts to deviate from that expected for an exponential 
disc. The error introduced by the subjectiveness of this pro- 
cedure is estimated by performing several 2D fits in which 
the fitting boundary is varied by 10-20 percent. 

3.3 Tests on artificial images 

In order to assess the reliability of the method and the effec- 
tiveness of the mask used at \z\ < 1.5 h z , a series of tests was 
performed on artificial images. First, we created fifty artifi- 
cial images using Eqs. (2) and (3). For the input parameters, 
random values were drawn from a uniform distribution cov- 
ering the following ranges: 
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Figure 1. The recovery of the disc parameters in tests of the 2D fitting method on artificial images (sec text), (a) — Edge-on central 
surface brightness, (b) - scalclcngth and (c) - scalchcight. 



Disc Bulge 

16 < no < 21 mag arcsec -2 17 < fj, c < 22 mag arcsec -2 

10 < hn < 40 arcsec 0.2 < r e < 10 arcsec 

2 < h z < 10 arcsec 0.5 < q < 1.0 

where /j,o an d /i e are used to denote Eo and S e in logarithmic 
units. To ensure the absence of artificial galaxies resembling 
bulge-dominated systems, the bulge-to-disc luminosity ratio 
was required to be less than two. Then, in order for the im- 
ages to resemble the /-band observations, a sky background 
corresponding to fi = 20 mag arcsec -2 was added. The use 
of a different value for the background in each image is un- 
necessary, because sufficient contrast with the background 
is already provided by the variation in fj,o and fie ■ The effect 
of seeing was mimicked by convolving the images with a cir- 
cular Gaussian of 1.5 arcsec FWHM, which is typical for the 
observations. Finally, Gaussian read-out noise and Poisson 
noise were added and the background was subtracted. 

For each artificial galaxy, initial estimates of the free 
parameters were generated by randomly offsetting the input 
parameters up to ±25 percent. The region corresponding to 
\z\ < 1.5 h z , was masked out. Then, the 2D fits were per- 
formed, while discarding data outside the la contour of the 
model and redefining the mask and the pixel weights in each 
iteration. Generally, convergence was achieved within three 
iterations. In many of these 2D fits the bulge parameters 
did not converge. This is mainly due the use of the mask, 
causing the bulge contribution in the fitted region to be too 
small (less than 10 percent). However, for the same reason, 
the 2D fits could be safely repeated by fitting a disc-only 
model (cf. Fig. 1). 

Figure 1 displays the test results for the disc parame- 
ters. It shows that the scale parameters are always recovered 
to within 10 percent, the standard deviation being about 2 
percent for each parameter. As expected, the accuracy of the 
recovery decreases with decreasing central surface bright- 
ness. For example, the three outliers of Fig. lb are all discs 
of low central surface brightness, fxo > 20 mag arcsec -2 . 
The amplitude of the deviations shown in Fig. 1 is mainly 
caused by the application of the mask; if no mask is applied 
the errors are less than 1 percent. Note that the application 
of the mask does not introduce any systematic error. 



We repeated the tests after incorporating a strong dust 
lane. For the distribution of the dust, we chose a double ex- 
ponential With /lR,dust = hn, stars and h Zjdust = 0.5 ft z , stars 

neglecting scattering (see Kylafis & Bahcall 1987, their 
cq. (14)). Of course, in reality, dust shows clumpy and fil- 
amentary structures on small scales. However, the global 
nature of the fitting routine and the edge-on view largely 
average these out (Kuchinski et al. 1998), so that the ex- 
ponential distribution is a good global approximation. Our 
neglect of scattering means that the effect of dust is slightly 
overestimated, since scattering would reduce the extinction. 
To resemble a strong dustlane, we chose a face-on central 
optical depth tj(0) = 1.0 (cf. Xilouris et al. 1999). 

The effect of the dust on the scale parameters turned 
out to be moderate; the hn and h z were found to be on 
average 4 percent and 9 percent larger. The use of a smaller 
mask than that adopted here (or no mask at all) would mean 
that the effect becomes substantial, causing errors on the 
scale parameters larger than 10 percent. 

The effect of a moderate stellar warp was investigated 
by simply introducing a vertical offset z W arp (R, <j>) to the 
axisymmetric disc model (Eqn. 1) for R > i? W ar P : 

L{R,Z,<t>) = L(R) e -|*-*warp(fl,«|/h* ; 
Zwarp(-R, <t>) = W [(R - -Rwarp)//lR.] COs(^) ft z , (5) 

where w is the warping rate in units of scaleheight per scale- 
length, .Rwarp designates the onset of the warp, and the line 
of nodes of the warp is assumed to be straight and along the 
line of sight (such that the effect of the warp is the largest) . 
For any reasonable warping, the effect on the derived scale 
parameters turned out to be small. For example, for a stellar 
warp with -R W ar P = 3 hn and w = 0.5, hn and h z are respec- 
tively under-and over-estimated by only 2 percent. This is 
mainly due to the fact that the inner parts of the edge-on 
disc are mostly unaffected by the warping. 

We also studied the errors introduced by the fixed pa- 
rameters. This was done by offsetting these parameters by 
their typical uncertainty and repeating the 2D fits. The fol- 
lowing fixed parameters were investigated: sky subtraction - 
1 percent error, seeing - 0.l"error, central position - l"error, 
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Table 1. Results of the 2D least-squares fits in the / band for the disc component 

Columns: (1) Name (ESO-LV catalogue); (2) Hubble type, taken from the Lyon/Mcudon Extragalactic Database (LEDA). These types 
were assigned as outlined in de Vaucouleurs et al. (1991); (3) Heliocentric velocity corrected for Virgo-centric flow (LEDA); (4) /-band 
absolute magnitude (de Grijs 1998), calculated using column (3) and corrected for Galactic extinction; (5) Maximum rotational velocity 
from LEDA, obtained from a global Hi profile and corrected for turbulent motions. For the four galaxies, for which u ma x is not directly 
available (marked by an asterisk), its value was calculated from a least squares fit to the /-band Tully-Fisher relation; (6) Neutral hydrogen 
mass, calculated using the Hi flux density from LEDA and column (3); (7) and (8) Edge-on central surface brightness and error; (9) and 
(10) Scalclength and error, in arcsec; (11) Scalclcngth in kpc; (12) and (13) Scaleheight and error, in arcsec; (14) Scaleheight in kpc; (15) 
and (16) The ratio of scalclcngth to scaleheight, i.e. the disc flattening, and error 



Galaxy type v vil A/J "max A4hi ^q,i ± fog, ± hn h z ± h z h R /h z ± 

(km s _1 ) (mag) (km s _1 ) (10 9 Mq) (mag arcsec -2 ) (arcsec) (kpc) (arcsec) (kpc) 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) (16) 
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position angle - 1 degree error, and inclination - 3 degrees 
error. For example, for a disc with hn/h z = 7 and an in- 
clination of 87 degrees, our 2D method overestimates the 
scaleheight by 4 percent. The scalelength was found to be 
most sensitive to an error in the sky background, resulting 
in an error < 4 percent. We estimate the total error due all 
systematic effects at < 6 percent. 

Further tests confirmed that the choice for the model 
used for the vertical light distribution is rather arbitrary 
when the region corresponding to \z\ < 1.5 h z is discarded 
(as stated in Sect. 3.1). For example, when instead a vertical 
sech(z/h z ) distribution is used to create the artificial images, 
the scale parameters found by our 2D method differ from the 
input values by about 2 percent. 



The algorithm does not take into account the effect of 
seeing on the disc component. This omission was tested by 
creating a number of artificial images, each with a different 
PSF, and applying the 2D fitting method. For the unlikely 
case that the seeing FWHM is comparable to the scaleheight 
of the disc, the error made by applying the method is only 3 
percent. Since the scaleheight is much larger than the seeing 
FWHM for all galaxies in our sample, one can safely ignore 
this effect. 

Finally, we investigated the effect of a truncation by 
applying the method to truncated artificial galaxies. These 
were constructed by including, in Eq. (1), L(R,z) = 0.0 
for R > i?max and numerically integrating along the line of 
sight. The difference in edge-on surface brightness between 
such a truncated disc and the infinite exponential disc in- 
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creases with projected radius, causing the 2D method to un- 
derestimate the exponential scalelength. For galaxies with 
-Rmax/^R > 4 this difference is small. For example, in a 
galaxy for which i? max /ftR = 4, our method underestimates 
the scalelength by ~ 6 percent. This error increases slightly 
for galaxies with smaller i? max //iR, reaching ~ 9 percent for 

Rmax/ftR = 3. 

In summary, we conclude that the uncertainty in the 
scale parameters obtained with the 2D method is caused 
by several effects. In order of decreasing importance, these 
are: (1) the presence of a truncation, the error depending on 
the value of -R m ax//iR, (2) residual dust extinction (error < 9 
percent), (3) inaccurate sky subtraction (error < 4 percent), 
and (4) the application of the mask to \z\ < 1.5 h z (error 
~ 2 percent, see Fig. 1). For example, for discs with an 
Rmax/hn of 4 (Sect. 5.2), the total measurement error in 
the scalelength is about 10 percent. 



Table 2. Results of the 2D least-squares fits in the I band 
for the bulge component 

Columns: (1) Name (ESO-LV catalogue); (2) Bulge type: e = 
exponential Eq. (3), v = de Vaucouleurs Eq. (4); (3) and (4) 
Bulge effective surface brightness and error; (5) and (6) Bulge 
effective radius and error; (7) and (8) Bulge axial ratio and error. 
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4 RESULTS 

4.1 The two dimensional fits 

First, to prepare the images for the 2D fits, Hn regions, 
foreground stars and background objects were masked out 
using a combination of the SEXTRACTOR algorithm (Bertin 
& Arnouts 1996) and manual editing in the galaxy's imme- 
diate vicinity. The centre and major axis position angle of 
the galaxies were determined by fitting ellipses to the /-band 
images using the GALPHOT package in iraf (Franx, Illing- 
worth, & Heckman 1989) (see Table 6, appendix). Then, 
several profiles were extracted in the / band, both paral- 
lel and perpendicular to the major axis. The latter profiles 
were inspected for the effects of dust extinction. In most of 
these profiles, the turnover due to dust extinction lies below 
z — 1.0 h z (see appendix, and de Grijs et al. 1997). However, 
in 16 galaxies the extinction clearly extends further out on 
one side of the galactic plane, e.g. in ESO 435-G14. This 
indicates that these galaxies are not exactly edge-on, and in 
these cases the affected side was excluded from the 2D fit. 

Fits were made to each /-band image, taking the values 
determined in the ID analysis (de Grijs 1998) as initial es- 
timates. In 10 cases, the mask needed to be decreased from 
| z | < 1.5 h z to \z\ < 1.0 h z in order to achieve conver- 
gence (Table 6, appendix). A bulge-disc decomposition was 
attempted in the 23 systems which clearly show an addi- 
tional central component. As the best-fitting bulge model, 
either exponential (Eq. 3) or r 1 / 4 (Eq. 4), the bulge-disc fit 
with the smallest reduced \ 2 was retained. However, in 14 
galaxies the bulge parameters did not converge. This was 
not unexpected, since in many galaxies only a small por- 
tion of the bulge light is emitted in the region outside the 
mask. Still, two of the 14 systems do show an extended cen- 
tral component; ESO 321-G10 and ESO 416-G25. In these 
galaxies this component has a boxy or peanut-like shape, 
which probably caused the fit to diverge. To determine the 
disc parameters in these 14 cases, we assigned an inner radial 
fitting boundary, as in the 'marking the disc' method (Free- 
man 1970), and repeated the 2D fits with a disc component 
only. Tables 1 and 2 show the resulting parameters for the 
entire sample. The errors in the parameters are not the for- 
mal errors of the least-squares fit; these were in general less 



than a few percent and do not take into account the subjec- 
tiveness of the assignment of the radial fitting boundaries. 
Instead, the errors were estimated by comparing results from 
several similar fits in which the inner and outer radial fitting 
boundaries were varied by 10-20 percent. For each galaxy, 
the 2D fits are further illustrated in the appendix. 

The 2D fits were repeated in the B band, using exactly 
the same fitting regions to further investigate the effect of 
dust extinction. For this, we calculated the scalelength ra- 
tios, /ir/^r- A scalelength ratio is an indirect measure of the 
radial colour gradient of the disc. Since dust is known to be 
concentrated to the galactic centres and dust extinction in- 
creases towards bluer wavelengths, the presence of dust will 
result in a larger h^/h^. Consequently, when the observed 
scalelength ratio is compared to the ratio expected solely 
from radial stellar population changes, the effect of dust ex- 
tinction can be inferred (Peletier et al. 1994; de Grijs 1998). 
The value of the scalelength ratio due to population changes 
is, however, somewhat uncertain. Peletier et al. (1994) ar- 
gue that the scalelength ratio between the B and the / band 
due to population changes is in the range 1.05-1.1. From 
modelling of 6 nearby edge-on galaxies, which includes a ra- 
diative transfer description, Xilouris et al. (1997, 1999) find 
an average ratio of 1.25 ± 0.13 (lcr). Recently, from a study 
of 21 Sb galaxies of various inclination, Cunow (1999) de- 
rives a ratio of 1.13. For our sample, the average h^/h^ is 
1.17 ± 0.18 (lcr). This is clearly within the range of ratios 
attributed to radial population changes, and smaller than 
the ratios observed in less-inclined galaxies (Peletier et al. 
1994). We conclude that the dust extinction in the fitted 
regions is small. This is in accordance with the results of 
Sect. 3.3. 

In Fig. 2, the disc scale parameters are compared to 
those obtained from the ID method (de Grijs 1998). For the 
scalelength, the differences are significant. The scalelengths 
determined with the 2D method are on average 13 percent 
larger. These differences arise mainly due to the combination 
of three effects. First, the fitted regions are different in both 
methods. According to Knapen & van der Kruit (1991), a 
difference of about 10 percent is not unusual in this respect. 
For example, the ID method adopts a fixed radial outer 
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Figure 2. Comparison of the disc scale parameters obtained with 
the 2D (this paper) and ID (dc Grijs 1998) method, (a) - The 
relative difference between the scalclcngths as a function of the 
scalelength determined with the 2D method, (b) - as (a) but for 
the exponential scaleheight. A few galaxies for which there is a 
large discrepancy are indicated. 



fitting boundary at 4 scalelengths. However, many galaxies 
show a steepening of the radial surface brightness in this 
region, often caused by a truncation (Sect. 4.2). This effect 
causes the scalelength to be underestimated. Secondly, the 
bulge contribution is substantial in early-type spirals, even 
in the outer parts of the disc, also causing an underesti- 
mation of the scalelength. Thirdly, the ID method adopts 
an exponential as a fitting model, instead of its projection 
(Eq. (2)). This effect leads to an overestimate of the scale- 
length by about 20 percent (Fig. 3). 

For the vertical scaleheight, the differences between the 
two methods are mostly within 25 percent (Fig. 2b), ex- 
cept for 7 outliers for which the 2D method finds a scale- 
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Figure 3. The systematic error made in the ID method (dc 
Grijs 1998) by fitting an exponential to the data instead of 
Eq. (2), which leads to scalelengths being determined 20 percent 
too large, (a) — The radial surface brightness profile according to 
Eq. (2) (solid line). The exponential profile is shown for compari- 
son (dashed line, arbitrary offset), (b) - The ratio of the measured 
scalelength, obtained in an exponential fit to the region between 
R = R m i n and R = 4 h R (de Grijs 1998), to the actual scalelength 
as a function of R mln _. The vertical dotted line indicates R m i n = 
h R used by de Grijs (1998). 



height which is smaller by about 50 percent. These differ- 
ences are caused by the differences in the fitting methods; 
the 2D method uses a two-dimensional fitting region and 
a bulge-disc decomposition, whereas in the ID method the 
scaleheight is determined from profiles taken parallel to the 
major axis at various projected radii. 

Our sample has 5 galaxies in common with the sample 
of Pohlen et al. (2000b). The model used by these authors is 
similar to that used here, but the method is different at two 
important points; (1) - no mask is applied to the region near 
the galactic plane and (2) - whereas with our method the 
truncation is largely avoided, their method is based upon 
the assumption that all discs show an infinitely sharp trun- 
cation. Table 3 lists both, the /-band scale parameters and 
truncation radii (Sect. 4.2), and the same parameters de- 
rived by Pohlen et al. (2000b) in the r band. The scale 
parameters agree well, except for the scalelengths of two 
galaxies, ESO 269-G15 and ESO 564-G27. For these, the 
determination of Pohlen et al. (2000b) is larger by about 30 
percent, much larger than the difference expected between 
the I and the r-band from radial population changes alone. 
We note that for ESO 269-G15, our 2D fit was performed 
to the eastern side because it is not perfectly edge-on (see 
appendix). In ESO 564-G27, additional extended emission is 
present at large distances from the plane in the inner parts, 
indicating that the use of different fitting regions may be 
the cause of the discrepancy. 
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Table 3. Comparison of the disc parameters with the lit- 
erature for five galaxies 

Columns: (1) Name (ESO-LV); (2) Band, I = /-band result ob- 
tained in this paper, r = Thuan-Gunn r-band result obtained by 
Pohlcn ct al. (2000b); (3) and (4) Scalelength and error; (5) and 
(6) Scalcheight and error; (7) and (8) Truncation radius and error 
(Sect. 4.2). 
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4.2 Fitting the truncation 

The truncation, or cut-off, of a stellar disc can be inferred 
from a sharp drop in the radial light profile (e.g. KS1-3). 
The truncation radius is then obtained by extrapolating the 
radial light profile to zero flux. The error made in such an 
extrapolation is largely due to Poisson statistics, imperfect 
sky subtraction and foreground stars (KS1-3). While trun- 
cations are most easily observed in edge-on systems because 
of their higher surface brightness, there one views the result 
of a complicated line of sight integration across the trun- 
cated disc, which includes the effects of dust extinction and 
non-axisymmetric components such as spiral arms. Addi- 
tionally, the truncation may not be axisymmetric, such that 
in an edge-on spiral the truncation becomes apparent at a 
different projected radius on either side (de Grijs et al. 2001). 
These considerations, as well as the variety of shapes of the 
radial light profiles seen in truncated edge-on spirals (e.g. 
Naslund & Jorsater 1997; de Grijs et al. 2001), illustrate 
the difficulties faced in the construction of a model for the 
shape of the truncation from studies of edge-on spirals. For 
these reasons we did not include a truncation in the model 
adopted for the 2D fits (Sect. 3). 

By applying an edge-fitting method similar to the one 
introduced by de Grijs et al. (2001), it is possible to objec- 
tively determine the truncation radius. The adopted method 
is as follows. A radial profile is constructed by first averag- 
ing vertically over the range 1.0 h z < z < 2.0 h z on the 
least obscured side of the galactic plane. The choice for this 
range minimizes the effects of dust and young populations, 
while retaining a good signal-to-noise ratio. To further in- 
crease the signal-to-noise in the outer parts the profile is 
also averaged radially, by using a semi-logarithmic scheme 
(de Grijs et al. 2001). We made sure that this radial averag- 
ing does not smooth out the truncation by testing the edge- 
fitting method on artificial images (these were constructed 
as in Sect. 3.3 and included the infinitely sharp truncation). 
Then, the part of the profile containing flux from the galaxy 
is defined by clipping the vertically averaged profile at the 
la level (a being the standard deviation of the background 
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Figure 4. The edge-fitting method, illustrated with a close-up 
of the averaged profile on one side of an artificial galaxy. This 
artificial galaxy was constructed as in Sect. (3.3) with i? max = 
100 arcsec and Rmax/^iR = 4. (a) — Linear surface brightness 
as a function of radius after vertical and radial averaging (open 
circles). The profile expected for an untruncatcd exponential disc 
(Eq.(2)) is shown by the dashed line. The vertical dotted line 
marks the outermost radius for which data is included in the trial 
fits. In this example, the minimum reduced \ 2 is reached when 
including 10 points in the fit (after trying 4 points and 7 points); 
the corresponding best-fitting straight line is shown as a solid line. 
Extrapolation of this fit to zero flux yields R = 100.3 ±3.6 arcsec. 
The difference in surface brightness with the untruncated disc at 
the outermost radius is marked by the arrow. It corresponds to 
3a of the background, (b) - as (a) but on a logarithmic scale. In 
addition, the vertically averaged profile is shown, shifted by 1.5 
mag arcsec -2 . The horizontal dotted line, shifted by the same 
amount, marks the standard deviation of the background in the 
vertically averaged profile (see text for further details). 



in the vertically averaged profile), and using it as a mask 
(i.e. a method of conditional transfer). In this way, the out- 
ermost projected radius included in the profile is determined 
directly by the data quality. Following this, a series of trial 
fits of a straight line is performed in linear surface brightness 
at the edge of the profile (note that de Grijs et al. (2001) 
used an exponential function). In the first trial fit the four 
outermost points are included. Then, in each of the subse- 
quent trial fits the inner fitting boundary is decreased by 
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three points. Finally, only the fit with the smallest reduced 
X 2 is retained and extrapolated to zero flux. The correspond- 
ing projected radius is then the adopted truncation radius, 
provided it satisfies the criterion outlined in the following. 

An edge-fitting method was already applied, although 
in a somewhat different form, to four sample galaxies in B, 
V, and / band (de Grijs et al. 2001). There, as in all pre- 
vious studies, the presence of a truncation is deduced from 
a qualitative comparison between the observed radial pro- 
file and the exponential disc model. In our sample, which 
is expected to contain many less obvious cases, a quantita- 
tive criterion is needed to determine its presence. To search 
for such a criterion, we performed tests of the edge-fitting 
method on artificial images. Not surprisingly, it was found 
that -R m ax is more easily recovered when the difference be- 
tween the observed surface brightness and that expected for 
an untruncated exponential disc (Eq. (2)) is large. In partic- 
ular, when the difference between the observed and expected 
surface brightness at the outermost radius of the profile is 
larger than about 2a, then i? max could always be recovered 
(a is the standard deviation of the background in the orig- 
inal image). When the difference at the outermost radius 
is smaller than 2a, which is caused by either a low central 
surface brightness and/or large 7? ma x//iR, then the recovery 
often fails. In what follows, a profile is said to be truncated 
if it satisfies this 2a criterion. The edge-fitting method is 
further illustrated in Fig. 4. 

Of course, when measuring surface brightnesses in the 
outer parts of galactic discs, the accuracy of the back- 
ground subtraction becomes very important. The main fac- 
tors which limit this accuracy are large scale flatfielding er- 
rors and the presence of background sources and foreground 
stars; our small-scale flatfielding errors are within 0.5 per- 
cent (de Grijs 1997) and the CCDs were not affected by 
fringing. The background emission in the /-band images was 
originally estimated by fitting a plane to regions far away 
from the galaxy (de Grijs 1997). To check the influence of 
residual light due to objects within that region, we repeated 
these fits after heavily growing of our object masks. The 
results show that the earlier fits slightly overestimated the 
background, on average by 0.3 times the standard deviation 
in the background. Although this oversubtraction is small, 
much smaller than the amplitude necessary to produce ar- 
tificial truncations (see below), it would affect the deduced 
values for i? max and we have therefore decided to apply the 
edge-fitting method to the new background-subtracted im- 
ages. To estimate the amplitude of remaining variations due 
to large-scale flatfielding errors, we inspected the distribu- 
tions of pixel values in 10x10 sized boxes at a large number 
of positions in the sky-subtracted background images. The 
standard deviation of the medians of these distributions was 
adopted as the error in the background determination. 

The results of the application of the edge-fitting method 
to our sample are shown in Table 4 (and the appendix). To 
calculate the error on 7i max two additional images were con- 
structed for each galaxy by adding/subtracting the error in 
the background determination. The fits were repeated and 
the error on /? max was calculated by quadratically adding 
half the difference between the truncation radii found from 
the two modified images to the formal error of the origi- 
nal fit. For the galaxies for which no truncation could be 
determined, Table 5 shows a lower limit for /? max . 



Table 4. Results of the successful truncation fits in the / 
band 

Columns: (1) Name (ESO-LV), an asterisk indicates a galaxy 
which is not exactly edge-on and shows signs of spiral structure; 
(2) Side; (3) and (4) Truncation radius and error, in arcsec; (5) 
and (6) Truncation radius and error, in /-band scalclcngths; (7) 
Truncation radius in kpc; (8) Difference between the observed 
surface brightness and that of the untruncated exponential disc 
(Table 1) at the outermost point of the profile, in units of the 
standard deviation of the background. 
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N 


126 


9 


3.4 


0.3 


14 


13 


157-G18 


NE 


127 


18 


3.5 


0.5 


9 


5 




SW 


143 


18 


3.9 


0.5 


10 


4 


201-G22 


NE 


102 


5 


4.2 


0.5 


25 


4 




SW 


96 


7 


3.9 


0.5 


24 


6 


202-G35* 


NW 


95 


6 


4.1 


0.5 


10 


6 


240-G11* 


SE 


210 


14 


4.3 


0.8 


36 


11 


269-G15* 


N 


97 


18 


3.4 


0.9 


20 


5 


288-G25* 


NE 


91 


19 


5.0 


1.5 


14 


7 


315-G20 


NE 


71 


4 


3.1 


0.7 


21 


6 


321-G10 


NE 


65 


6 


3.2 


0.7 


12 


5 


416-G25 


NE 


74 


4 


3.3 


0.4 


23 


5 




SW 


74 


3 


3.3 


0.3 


23 


5 


435-G14* 


SW 


79 


12 


4.5 


1.1 


13 


6 


446-G18* 


s 


84 


12 


3.3 


0.6 


25 


13 


446-G44 


W 


79 


6 
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0.3 


14 


26 




E 


88 
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2.8 


0.3 


15 


25 


460-G31* 
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113 


24 


4.0 


1.0 


42 


7 


487-G02 


NE 
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11 


4.6 


0.6 


11 
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SW 


110 
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4.5 


0.4 


11 


10 


509-G19* 


SW 


90 


12 


3.6 


0.8 


62 


8 




NE 


91 


9 


3.7 


0.7 


62 


9 


564-G27* 


N 


148 


11 


3.9 


0.6 


19 


8 



To summarize, a truncation is found in 28 of the 68 pro- 
files (counting each side of each galaxy separately). At the 
outermost point in these profiles, the differences between the 
exponential disc model and the data are larger than ~ 3a, 
clearly satisfying our criterion. These differences are also 
much larger than the error in the background determination, 
i.e. the truncations are not produced by over-subtracting the 
sky background (Table 4 and the appendix). Of the other 
40 profiles, about half had to be discarded because of the 
presence of foreground stars, background objects or a lim- 
iting field of view. In the remaining profiles, the difference 
between the data and the untruncated exponential disc at 
the outermost point is less than 2a. The 28 truncations are 
found in 20 galaxies, 8 of which show a truncation on both 
sides. In all of these 8 cases the truncation appears sym- 
metric; any asymmetries are less than 15 percent. Of the 12 
galaxies in which a single truncation is seen, there are only 2 
for which the profile at the other side of the galaxy could be 
studied as well; ESO 435-G14 and ESO 460-G31. Looking 
directly at the images of these, we note that their discs are 
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Table 5. Lower limits to the truncation radii for the un- 
successful truncation fits 

Columns: (1) Name (ESO-LV); (2), (3) and (4) Lower limit to the 
truncation radius in arcsec, /-band scalelengths and kpc respec- 
tively 



Oalaxy 


(-Rmax)mi 






(arcsec) 


(%) 


(kpc) 


(i) 


(2) 


(3) 


(4) 


033-G22 


80 


4.4 


21 


138-G14 


130 


3.1 


13 


263-G15 


113 


3.6 


17 


263-G18 


80 


2.9 


19 


322-G87 


60 


3.6 


13 


340-G08 


70 


6.1 


12 


340-G09 


70 


3.1 


11 


435-G25 


200 


2.1 


30 


435-G50 


60 


4.5 


10 


437-G62 


115 


3.2 


21 


506-G02 


65 


3.4 


16 


531-G22 


85 


4.3 


18 


555-G36 


65 


5.0 




575-G61 


55 


3.4 


6 



not perfectly edge-on and show a signature of strong spiral 
arms which may be causing the single-sided truncation (cf. 
de Grijs et al. 2001). Galaxies which are not exactly edge- 
on and show signs of spiral structure are marked with an 
asterisk in Table 4. As an alternative possibility, the single 
truncations may be caused by a non-axisymmetric or lop- 
sided stellar disc in which the truncation occurs at a smaller 
projected radius (= higher surface brightness) on one side. 

To assess the effect of dust extinction, the fits were re- 
peated using the region below one scaleheight. The resulting 
Umax are consistent with the values of Table 4, indicating 
that dust extinction near the edge of the disc is not very 
important. Observations of dust tracers such as CO emis- 
sion and mm/submm continuum in the outer parts of spiral 
galaxies confirm this picture (Combes & Becquaert 1997; 
Neininger et al. 1996). 

Other effects which may cause an apparent truncation 
are warping or flaring of the stellar disc. Although our sam- 
ple contains no spirals which are strongly warped (Sect. 2), 
artificial images of warped discs (Sect. 3.3) show that a mod- 
erate warping can still introduce a break in the radial profile 
similar to the signature of a truncation. More importantly, 
away from the major axis the radial profiles show a strong 
asymmetry; a break at a small galactocentric radius on one 
side and an upturn followed by a break at a large galacto- 
centric radius on the other. This causes a large difference 
between the apparent truncation radii of both sides. While 
similar upturns in a few of our sample galaxies may be a 
sign of warping (or the effect of spiral arms), the galaxies 
having a double-sided truncation show no evidence for such 
a strong asymmetry. To address disc flaring we also investi- 
gated artificial images of flared discs viewed edge-on. These 
which were constructed as in Sect. 3.3, but included a scale- 
height increasing linearly with radius. It turned out that a 
very strong and localized thickening in the outer parts is re- 
quired to produce a feature resembling a radial truncation 
(a rate larger than 1 h z per hn). For our sample galaxies, 
and late-type spirals in general, such thickening is not ob- 



served (de Grijs & Peletier 1997; Morrison et al. 1994; Fry 
et al. 1999). We conclude that warping and flaring are very 
unlikely to have caused the truncations. 

Again, we can compare our results to those obtained by 
Pohlen et al. (2000b), who estimated the truncation radius 
by visually comparing radial profiles to a truncated model 
(Table 3). For these 5 galaxies the results for the trunca- 
tion radius are consistent. Thus, the truncation is seen in 
different data-sets using different methods, indicating that 
the truncations are real. This is further emphasized by the 
tight relation observed between hn and i? ma x (Sect. 5.2). 

Finally, we would like to emphasize that the fraction 
of galaxies in which a truncation is found, about 60 per- 
cent, is a lower limit. Clearly, the detection of a truncation 
depends on the edge-on central surface brightness, on the 
Rmax/h-R of the disc, on the limiting magnitude reached by 
the observation and, to a lesser extent, on the shape of the 
truncation. Since these elements are different in each image, 
the selection effect is different in each case. 



5 DISCUSSION 

5.1 The disc flattening 

We further investigate the disc scale parameters obtained 
in Sect. 4.1. In Figs. 5a &d the distributions of the /-band 
scalelengths and scaleheights are plotted. The corrected dis- 
tributions were obtained by applying a volume correction 
to the diameter-limited distributions (van der Kruit 1987; 
Davies 1990; de Jong 1996b) and normalizing to the num- 
ber of galaxies in the sample. The diameter selection limit 
implies that at small physical sizes the sampled volumes are 
too small to contain any edge-on spirals (from our sample 
the density of spirals with inclination > 87 ° is only about 2 
x 10 Mpc~' ). This selection effect causes the low number 
of small spirals seen in the scalelength distribution. Small 
scalelength spirals are known to dominate the population of 
spirals (van der Kruit 1987; de Jong 1996b) . Recognizing the 
fact that these small spirals are underrepresented, the shape 
of the volume-corrected distribution is consistent with the 
results of van der Kruit (1987) and de Jong (1996b). Interest- 
ingly, the distribution of scaleheight shows a similar shape, 
including a lack of small scaleheight galaxies (Fig. 5d). This 
lack is also due to the selection against intrinsically small 
galaxies because scaleheight and scalelength are observation- 
ally related (see below). Still, the smaller scaleheight discs 
are the most numerous; the volume-corrected distribution of 
scaleheights implies that about 90 percent of the spirals has 
h z < 0.6 kpc. 

For 33 spirals with known maximum rotational veloci- 
ties (Table 1), both scale parameters are plotted as a func- 
tion of maximum rotational velocity in Fig. 5b &e. This 
shows clearly that the disc scalelength increases with max- 
imum rotational velocity. A Spearman rank correlation test 
results in a coefficient of 0.75, corresponding to a confidence 
level greater than 99 percent. Since the maximum rotational 
velocity is an indicator of total mass, this merely shows that 
more massive spirals tend to have larger stellar discs. More 
interestingly, the disc scaleheight shows a similar behaviour 
(Spearman correlation coefficient 0.74). Hence, more mas- 
sive spiral galaxies have both larger and thicker discs. The 
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Figure 5. The disc scale parameters in physical units; (a) - The scalclcngth distribution before (dotted histogram) and after (solid 
histogram) the volume correction (see text), (b) - Scalclcngth as a function of maximum rotational velocity, (c) - Scalclcngth as a 
function of total Hi mass, (d)&(f) — As (a)&(c) but for the scaleheight. The dashed line in panel (e) is the relation expected for a 
vertically exponential disc with face-on central surface brightness fis = 21.65 mag arcsec~ 2 and constant mass-to-light ratio M/Lb = 2 
(M/L)q which obeys the vertical Jeans equation and Eq. (6). In panels (b), (c), (c), and (f) the pentagrams indicate the galaxies from 
KS1-3. 



most extreme example of this is ESO 509-G19. This huge 
spiral has an extremely large maximum rotational veloc- 
ity of 392 km s _1 and a disc scalelength and scaleheight 
of 17 kpc and 1.7 kpc respectively. A similar correlation 
is seen between the scale parameters and the total H I mass 
(Fig. 5c & f), with the notable exception of ESO 437-G62 (= 
NGC 3390). This outlier has a remarkably small amount of 
H I and a thick disc, and is classified as an SO-a in the ESO- 
LV catalogue. For spirals however, more massive Hi layers 
tend to be associated with larger and thicker stellar discs. 
We observe no clear trend with morphological type. 

Observations of the stellar kinematics in 12 spirals by 
Bottema (1993) have revealed a related trend; more mas- 
sive spirals tend to have larger stellar velocity dispersions. 
A simple straight line fit, adopted from his fig. 6, gives: 

Cz|r=o = (0.29 ± 0.10) (6) 

where a z |r=o is the central vertical velocity dispersion of the 
disc stars. This vertical velocity dispersion is directly related 
to the scaleheight of the disc mass distribution through the 
vertical Jeans equation. When the vertical velocity disper- 
sion is integrated in the vertical direction (see van der Kruit 
1988), one has: 



o- z | R= o = VCGEohz, (7) 

where Eo is the central disc surface density and C is a 
constant which depends on the exact choice of the verti- 
cal density distribution; for an exponential distribution C = 
37r/2, for a sech- dependence C = 1.7051 it and for the locally 
isothermal distribution C = 1. If we assume that the scale- 
height of the disc mass distribution is proportional to that of 
the old stellar population, Eq. (7) predicts that an increase 
in the vertical velocity dispersion is associated with an in- 
crease in the scaleheight of the light of the old population. 
Furthermore, by eliminating a z \R=o between Eq. (6) and (7), 
the scaleheight of the disc light is then also expected to be re- 
lated to the maximum rotational velocity; h z oc u 2 aax /\/So. 
As an example, we illustrate this in Fig. 5e (dashed line) 
for vertically exponential discs with face-on central surface 
brightness fis = 21.65 mag arcsec -2 (Freeman 1970) and 
constant mass-to-light ratio M/Lb = 2 (M/L)q. While the 
disc mass-to-light ratios in and among spiral galaxies are 
still poorly known, at least the shape of the predicted rela- 
tion for a constant mass-to-light ratio is roughly similar to 
that of the observations. 

In Fig. 6 we investigate the ratio of the /-band scale 
parameters, which is a direct measure of the intrinsic flat- 
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Figure 6. The disc flattening; (a) — The flattening distribution before (dotted histogram) and after (solid histogram) volume correction, 
as described in the text, (b) - Scalelength versus scalchcight, (c) - Flattening as a function of maximum rotational velocity, (d) - 
Flattening as a function of total Hi mass. In panels (b), (c) and (d) the pentagrams indicate the galaxies from KS1-3. 



tening of the disc light. The volume corrected distribution of 
the disc flattening (Fig. 6a) is remarkably narrow; about 70 
percent of spirals seem to have a flattening between 6 and 8. 
The average flattening of the diameter-limited distribution is 
(tin/hz) — 8.5±2.9 (lo"), while that of the volume corrected 
distribution is (/ir,//i z ) = 7.3 ± 2.2 (la). Figure 6b reveals 
a clear correlation between /ir and h z (Spearman correla- 
tion coefficient 0.84). Considering the errors, the spread in 
the /iR-fez plane is real, suggesting that the physical relation 
between /ir and h z cannot be described by a simple linear 
functionality expected for discs of constant flattening. This 
is even more clear in Fig. 6c &d. The range and maximum 
of the disc flattening seem to increase with both maximum 
rotational velocity and total H I mass. A similar behaviour 
is seen with absolute magnitude M t ° (not shown). The disc 
flattening shows a similar, but weaker, trend with Hubble 
type (de Grijs 1998; Schwarzkopf & Dettmar 2000), in ac- 
cordance with the results on the disc axial ratio obtained by 
Guthrie (1992) and Kudrya et al. (1994). On the other hand, 
the minimum disc flattening seems to be independent of all 
of these parameters. At present, there is no clear physical 
explanation for these trends. 

The disc flattening can be used to estimate the contri- 



bution of the disc to the rotation curve. For a self-gravitating 
exponential thin disc, there is a relation between maximum 
rotational velocity (i) and disc scalelength (Freeman 
1970): 

£ = 0.88 VttGEo/ir. (8) 

When the central surface density is substituted from Eq. (7), 
a simple relation between the maximum rotation, vertical ve- 
locity dispersion and flattening of the disc is obtained (Bot- 
tema 1993): 

W = (0.69 ± 0.03) a z | R=() ,/^, (9) 

V n z 

where the constant accounts for density laws ranging from 
an exponential to a sech-dependence, and hn/h z is the flat- 
tening of the disc mass. This equation states that a flat- 
ter disc with the same vertical stellar velocity dispersion or 
'temperature' will be more massive. By eliminating <r z be- 
tween this equation and the empirical relation Eq. (6), one 
can estimate the contribution of the disc to the rotation 
curve: 

«™/w = (0.21 ± 0.08) S^. (10) 
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Figure 7. (a) - The truncation radius versus disc scalelength. The dotted line is the least-squares fit of a straight line. The arrows 
indicate a lower limit for 13 galaxies for which no -R max could be determined. The pentagrams indicate the galaxies from KS1— 3. (b) 
— Kmax/'iR as a function of disc scalelength. The arrows and pentagrams are as in (a). The lines illustrate the prediction of a collapse 
model for disc galaxy formation (see text). Each line connects models of the same total mass, but with a varying spin parameter, A, 
ranging from 0.01 (at the left end of each line) to 0.28 (at the right end of each line). The total mass increases in powers of ten, from 
10 10 Mq (for the line at the left) to 10 13 Mq (for the line at the right). 



By taking hn/h z = 10 for the average flattening of the 
disc mass and using a locally isothermal vertical distribu- 
tion, Bottema (1993) argued that the disc supplies on av- 
erage 63 percent of the observed maximum rotation, i.e. 
Wmax/vmax = 0.63. Of course, this argument can be used 
for our sample, where we measured the flattening of the disc 
light in the / band. Assuming that the flattening of the disc 
mass equals that of the disc light from the old stellar pop- 
ulation, and using the volume corrected average of the flat- 
tening of 7.3 ± 2.2, we find (v££/iw) = 0.57 ± 0.22. This 
confirms the result of Bottema (1993) and suggests that in 
a typical spiral galaxy, the actual disc mass-to-light ratio is 
about one third of that expected in a maximum disc situa- 
tion (van Albada et al. 1985; Freeman 1992). However, the 
result is not entirely unambiguous. Given the error in the 
determination, the value of the disc contribution to the ro- 
tation could be about 80 percent. This would mean the disc 
can still be close to maximum, since, to avoid hollow halo 
cores and allow for a bulge contribution, a maximum disc 
has Wmax/vmax ~ 0.9. The relatively large error in our de- 
termination stems from the scatter in the empirical relation 
Eq. (6), a relation which is still ill-constrained due to the 
small number of spirals for which stellar velocity dispersions 
are known. 

5.2 The disc truncation 

The truncation radius is an important global parameter, for 
which we need to find an explanation in theories of disc 
galaxy formation and/or evolution. In Sect. 4.2, we showed 
that the light distribution of at least 20 out of the 34 galax- 
ies in our sample is radially truncated. These galaxies dis- 
play a remarkably tight relation between the truncation ra- 
dius and disc scalelength (Fig. 7a) (Spearman correlation 
coefficient 0.95). By taking the average ii max of both sides 
in case of a double-sided truncation, we obtain an average 
Rmax/hn = 3.6 ± 0.6 (la) in the / band. Previous studies 
have also revealed the relation between the truncation radius 



and disc scalelength, but yielded different values for the ratio 
-Rmax//iR. For example, for 7 nearby spirals KS3 find a mean 
value -Rmax//iR = 4.2 ± 0.5, Pohlen, Dettmar, & Liitticke 
(2000a) find 2.9 ±0.7 from 31 galaxies, whereas Schwarzkopf 
& Dettmar (2000) obtain 3.7T0.9 from 61 galaxies (although 
the latter two samples have several galaxies in common). 
These differences may be partly due to the use of different 
passbands; the use of bluer passbands will result in larger 
scalelengths, and hence leads to smaller i? max /ftR. 

We further investigate the distribution of -Rmax/ZiR in 
Fig. 7b by plotting 7? max //iR as a function of scalelength. 
For those galaxies for which no truncation was found, a lower 
limit is shown (arrows). The plot seems to reveal an increase 
in -R m ax//iR towards small scalelengths; the average ratio for 
the spirals with /ir < 4 kpc is 4.4, more than one standard 
deviation higher than that of the entire sample. This in- 
crease of -Rmax/^R may be related to its decrease found at 
very large scalelengths by Pohlen et al. (2000a). However, 
the reality of this feature is not clear considering the modest 
sample size and the selection effect against galaxies of small 
physical sizes. In addition, the ability to find a truncation 
depends, among other things, on the edge-on central surface 
brightness (Sect. 4.2). We are, for example, still missing low 
surface brightness galaxies, which may have entirely differ- 
ent iimax//lR- 

While our view of the distribution of -R m ax//iR is cer- 
tainly obscured by these complex selection effects, additional 
information can be obtained from face-on samples. We have 
estimated a minimum i? max for the spirals in the face-on 
sample of de Jong & van der Kruit (1994). By taking the 
lowest contour in their i?-band contour plots, we find an av- 
erage (i?max/^R) min = 4.0 ± 1.1 (lcr). When combined with 
the results from our truncation analysis, this suggests that 
the ratio of truncation radius to disc scalelength in small 
scalelength spirals is at least four. This is important, be- 
cause small spirals are the most numerous in the local Uni- 
verse (van der Kruit 1987; de Jong 1996b, Sect. 5.1). Note 
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that this also implies that the truncation is more difficult to 
detect in small scalelength spirals. 

We can compare the results on the truncation radii to 
the predictions of the collapse model of disc galaxy forma- 
tion (White & Rees 1978; Fall & Efstathiou 1980; Dalcanton 
et al. 1997). The collapse theory is compelling because its 
predictions are in general agreement with the basic structure 
and rotation curves of disc galaxies, over the full range of 
central surface brightness (Dalcanton et al. 1997; Mo, Mao, 
& White 1998). If it is further assumed that proto galaxies 
are sharp-edged then the collapse theory also predicts the 
outermost radius of the baryonic disc, which corresponds to 
the material with the highest specific angular momentum. 
To quantify that prediction and make a comparison with 
the observed truncation radii of stellar discs, we calculated 
model surface density profiles over a range in total mass, 
Aftot = 10 10 -10 13 M , and spin parameter, A = 0.01-0.28, 
using the method of Dalcanton et al. (1997). In the calcula- 
tions, we used a fixed baryonic mass fraction F = 0.10, and 
a Hernquist halo at the onset of baryon-decoupling (as in 
Dalcanton et al. 1997). From the calculated surface density 
profiles, the scalelengths were derived using a method simi- 
lar to the 'marking-the-disc' method (Freeman 1970) and the 
outermost radii were obtained by taking the radius at which 
the density drops to zero. The result is plotted in Fig. 7b as 
lines of constant mass. Both the outermost radius and the 
scalelength of the baryonic proto-disc increase with the mass 
and angular momentum of the proto galaxy, such that their 
ratio remains roughly constant at 3-4. Taking a different 
or non-constant baryon mass fraction does not significantly 
change this result. The predicted ratio is somewhat smaller 
than the ratio of 4.5, predicted by van der Kruit (1987) based 
on a comparison of the angular momentum distribution of 
an exponential disc to that of a uniformly rotating, uniform 
sphere with A = 0.07. 

The predictions of the collapse model of disc galaxy for- 
mation roughly coincide with the observed truncation radii 
of stellar discs, although they do not show an increase of 
Rmax/hR toward small scalelengths. It is clear from Fig. 7b 
that some galaxies fall above or below the predicted range, 
which may be attributed perhaps to environmental effects 
such as accretion and mergers. At any rate, in general the 
ratio Rmax/hn predicted by the collapse model is similar 
to that seen in the stellar luminosity density. This may in- 
dicate that the present radial structure of the stellar disc 
still largely traces that of the baryonic proto-disc. If true, 
this result implies that Hi present beyond the stellar disc 
(Broeils & Rhee 1997) is of a more recent origin (van der 
Kruit 1987). The result does not seem to be in accordance 
with the radial colour gradients observed in nearby face- 
on spirals. These indicate a stellar population change, with 
both the average stellar age and metallicity decreasing to- 
wards larger radii (de Jong 1996c). However, whether this 
shows that spiral galaxies have grown since their formation 
is largely an open question, for it is unclear if the radial 
colour gradient applies to the entire stellar disc or is caused 
by a young population in the plane of the galaxy. 

It should be noted that our observations compare 
equally well to the value predicted by the theory of stochas- 
tic self-propagating star formation, for which, assuming a 
flat rotation curve, Rmax/hn = 4 (Seiden et al. 1984). The 
other possible scenario, in which the truncation of the stel- 



lar disc is caused by a threshold on star formation due to 
its dynamical stability (Fall & Efstathiou 1980; Kennicutt 
1989), does not make a direct prediction for 7? max . How- 
ever, in a disc galaxy formation model similar to that of 
Dalcanton et al. (1997), van den Bosch (2001) has incor- 
porated a star formation recipe which includes a thresh- 
old surface density based on Toomre's (1964) criterion for 
local stability. For models of small bulge-to-disc ratio and 
with maximum rotational velocities exceeding 100 km s _1 , 
comparable to the galaxies considered here, the presence of 
this threshold produces truncations in the stellar discs with 
i?max/^R = 2-4 (his fig. 7). This is in good agreement with 
the observations and also corresponds to the outermost ra- 
dius of the proto-disc according to the collapse model inves- 
tigated above. In particular, this latter coincidence seems to 
be caused by the fact that for massive late-type spirals the 
model of van den Bosch (2001) predicts that star formation 
occurs essentially throughout the entire gas disc, such that 
the present day truncation radius of the stellar disc is only 
marginally smaller than the outermost radius of the gaseous 
disc. We tentatively conclude that although the origin of the 
truncation of the stellar discs is still unclear, and may well 
be caused by a threshold on star formation, the truncation 
radius mainly reflects the material of the highest specific 
angular momentum which has collapsed during formation, 
at least for the intermediate to late-type galaxies considered 
here. To test whether the stellar truncation radius corre- 
sponds to the predicted threshold radius on star formation 
would require a detailed study of H I and H2 surface density 
profiles, as first suggested by Kennicutt (1989). 



6 CONCLUSIONS 

We have analysed /-band photometry of 34 edge-on spiral 
galaxies, taken from the sample of de Grijs (1998), to study 
the flattening and truncation of stellar discs. The exponen- 
tial scale parameters of the disc light were derived using a 
2D least-squares fitting method. This procedure includes a 
bulge-disc decomposition and adopts a mask in the region 
\z\ < 1.5 h z to systematically exclude the region affected by 
dust extinction and young populations. It is better suited 
to extract the disc scalelength than the earlier ID analysis 
(de Grijs 1998). We used an edge- fitting method, similar to 
the one introduced by de Grijs et al. (2001), to study the 
disc truncation. This method uses a simple straight line as a 
fitting function and a quantitative criterion to establish the 
presence of a truncation. 

We find a clear increase in the scaleheight of the stel- 
lar disc as a function of maximum rotational velocity and 
total Hi mass. This is in general accordance with observa- 
tions of the stellar kinematics in spiral galaxies (Bottema 
1993); larger discs are clearly more dynamically evolved and 
thicker. Both, the maximum and the range of the flatten- 
ing of stellar discs seem to increase with maximum rotation 
and total Hi mass. We used the average volume corrected 
flattening of the disc light, {hn/h z } — 7.3, to estimate the 
disc contribution to the rotation curve, with the assumption 
that the flattening of the disc mass equals that of the disc 
light. This resulted in a disc contribution of 57±22 percent, 
which still - barely - allows for a maximum disc situation. 

At least 20 galaxies in our sample (60 percent) are 
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truncated, displaying a tight relation between scalelength 
and truncation radius. The average truncation radius corre- 
sponds to 3.6 ± 0.6 /-band scalelengths. In addition, there 
seems to be a small increase in 7? m ax/^R towards galax- 
ies of small scalelength. Small scalelength spirals, which are 
the most numerous spirals in the local Universe, have an 
-Rmax/ZiR of at least four. A comparison of this ratio to that 
predicted by theories of disc galaxy formation seems to indi- 
cate that for the intermediate and late-type galaxies consid- 
ered here the truncation radius mainly reflects the material 
of the highest specific angular momentum which has col- 
lapsed during formation. These scenarios are inconclusive 
with respect to the question of the origin of the truncation, 
and imply that a structure analysis of the truncation alone 
is insufficient to solve this issue. 

It is clear that many questions are left unanswered. In 
an attempt to partly remedy this situation, we are currently 
studying the stellar and H I kinematics of the galaxies in our 
sample. This will, for instance, further improve the statistics 
on the empirical relation between maximum rotation and the 
stellar velocity dispersion, which could significantly advance 
our knowledge on the dark matter content of spiral galaxies. 
Still, studies of larger samples are needed to further con- 
strain many of the relations found here, e.g. the trend of the 
disc flattening with disc mass. These future studies should 
preferably use deep, wide field, near infrared photometry. 
By targeting face-on as well as edge-on spirals, it becomes 
possible to more accurately study the disc truncation as well 
as its shape. 
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APPENDIX 

We further illustrate the results of the 2D fits (Sect. 4.1) 
and the truncation analysis (Sect. 4.2). For each galaxy, in- 
formation is presented in three panels. 

The panel at the top left shows the surface bright- 
ness distribution in the / band, after editing and rotation 
(Sect. 4.1). The surface brightness at the lowest contour 
corresponds to the 3a noise level of the background. The 
3(J level is marked by the solid arrow on the magnitude 
scale of the lower left panel. The contour interval is 0.5 mag 
arcsec" 2 , except for ESO 033-G22, ESO 041-G09 and ESO 
138-G14, where we chose an interval of 1.0 mag arcsec" 2 
for clarity. The orientation of each image is given in Ta- 
ble 6, columns (5) & (6). The inner and outer radial fitting 
boundaries and the mask (Sect. 4.1) are marked by the ver- 
tical and horizontal dotted lines respectively. The dashed 
contour marks the boundary outside of which the data were 
not included in the 2D fit (Sect. 3.2). 

The bottom left panel shows, on the same scale as the 
top left image, two profiles taken parallel to the major axis. 
In each case the best-fitting model is shown for comparison 
(dashed lines) . The bottom profile (solid line) is the average 
of two profiles. These were taken at a distance of 1.5 h z on 
each side of the galactic plane and averaged over a range of 
1/3 ft z to yield an acceptable signal-to- noise. When one side 
of the galaxy was fitted, only that side is used to construct 
this profile. The inner and outer radial fitting boundaries 
are indicated by the horizontal dotted lines. The top profile 
in this panel is the one used in the truncation analysis. It is 
shifted upwards by 1.5 magnitudes for clarity. If a truncation 
is present then the truncation fit (solid line) is also shown. 
The dashed-dotted lines indicate an envelope corresponding 
to ±3 times the background error. The open arrow at the 
magnitude scale, also shifted by 1.5 mag, indicates the stan- 
dard deviation of the background in the vertically averaged 
profile used to clip the light profile (see Sect. 4.2 for further 
details). 

The panel at the bottom right shows two profiles taken 
perpendicular to the major axis. The corresponding pro- 
jected radii at which these were taken are mentioned in Ta- 
ble 6, columns (7) & (8). In each case, the model is shown for 
comparison (dashed lines) . Again, each of the profiles is an 
average of two profiles, one taken on each side of the galaxy 
after averaging over a range of 1/3 /ir. The boundary of the 
mask is indicated by the dotted lines. 



Table 6. Additional properties of the 2D fits presented 
in the appendix 

Columns: (1) Name (ESO-LV), the galaxies for which a mask 
\z\ < 1.0 h z was used are indicated by an asterisk; (2) and (3) 
Centre used in the 2D fits; (4) Adopted major axis position angle 
(N — > E) (5) and (6) Image orientation; (7) and (8) Radii used to 
extract the vertical profiles. 
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